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1. Introduction 

Bezier curves and tensor-product Bezier surface patches have been used 
extensively in the computer graphics/computer aided design community for the 
modeling of physical objects for automated design. In 1980, Lane and Riesenfeld 
[Lane-Riesenfeld, 1980] presented algorithms for the evaluation of Bezier curves and 
surfaces, and showed how these algorithms could be used in solving many of the 
geometric problems in CAD, and how they could be used in the computer graphic 
rendering problems for such surfaces. These algorithms were based upon an elegant 
subdivision technique, which was developed directly from the mathematical properties 
of the Bezier curve. This technique allows the manipulation of a curve/surface 
without the need to evaluate the blending functions and their derivatives. 

Triangular Bezier patches have also been studied for the same reasons. Farin 
[Farin, 1982] has used them successfully in the interpolation of scattered 3-d data and 
Petersen [Petersen, 1984] has recently exhibited their use in a contouring algorithm. 
Subdivision of these patches has centered on a method [Farin, 1982] which 
subdivides the patch into three subpatches using intermediate control points 
calculated by the de Casteljau algorithm [de Casteljau, 1959]. However, due to the 
fact that one side of each of the three subtriangles is one of the sides of the parent 
triangle, the sides of the resulting subtriangles do not reduce uniformly in length. 

This implies that the polygonal rendering algorithms [Lane-Carpenter, 1979; Blinn, 
Carpenter, Lane, Whitted, 1980] (which require surfaces to be subdivided until they 
are flat) cannot be used with this subdivision method. (It is noted that Petersen used 
this method in his contouring algorithm, however, an extra subdivision step was 
added to insure that the lengths of the sides reduced uniformly.) 

This paper presents a subdivision method for Bezier triangles that parallels the 
development of Lane-Riesenfeld for Bezier Curves. The method is developed from the 
mathematical description of the patch, and the resulting algorithm does not require 
evaluation of the blending functions. This method subdivides the triangle uniformly 
so that the lengths of the sides of the subtriangles are uniformly reduced. This 
method is a member of the class of methods first presented by Goldman [Goldman, 
1983], which contain subdivision methods that arise from degenerate subdivision of 
Bezier tetrahedra. In this paper, we show that a certain method is a natural 
extension of the Lane-Riesenfeld results for Bezier curves. 

2. Bezier Triangles 

The Bezier triangle of degree n , defined by the set of control points 
P = {Pi j k : t+/+A=n} over the triangle T is given by 
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where 


B„{P]T)= Yi uP2,Ps) 
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are the bivariate Bernstein polynomials of degree n , using the barycentric coordinate 
system defined on the triangle T . 


3. Subdivision of Bezier Triangles 

Consider the triangle T with vertices P 1 P 2 , P 3 , where P t t R 2 . Let T be 
parameterized by the barycentric coordinates Consider the subdivision of 

T given by the following sets 
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|p e T : P — Pi P i+P 2 P 2+M3 P 3 ) 
| P f T : P = p 1 P ^p^P 2 +p 3 P 3 , 
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| P < T : P = p x P x +p 2 P 2 +p 3 P 3 , 
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The subdivision algorithm then proceeds as follows. 


Subdivision Algorithm 

Let P = {P t j * : » +/ +/: =n } be a polygonal mesh defined over the triangle 
T and let 
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if n j > 0 


if n 2 > 0 


if n 3 > 0 

ifn 1 = n 2 = n 3 = 0 


Then 
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B„(P ; T) = B n ( Qi ; T t ), for = 1, 2, 3, 4 

where 

Q i = jp.-fo 0 ’* 1 :•+*=», m = 0, .} 

Q 2 = |^,!y’,o’ 0i : *+/=»> m =0, jj 

Q 3 ==z : ;'+*=«, m = 0, .... *} 

Qa = {^V* 1 :»+/+*=»} 

Lemma 

Let Q = |§ t ; ^ : * + /+A? = n| be the control points for the subdivided 
Bezier triangle over T 4 ; then we have 

I Qi+l,j ,k ~Qi,j,k I — M 

1 Qi,j+l,k~Qij ,k I — M 
I Qi,j ,k+l~Qi,j ,k I — M 

where 

j I P i+i ,j,k-P » ,j,k I 

M = —max | Pi j +i i ~Pi j k I : 1 + / = n ,andt ,j,k <n 

2 [ | Pi] i 9 k+i-Pij'.t I 


69 



4. Bibliography 


[Bezier, 1972]. Bezier, P.E., Numerical Control -- Mathematics and Applications, John 
Wiley and Sons, London, 1972. 

[Blinn, Carpenter, Lane, Whitted, 1980]. Blinn, J.F., L.C. Carpenter, J.M. Lane, and T. 
Whitted, Scan Line Methods for Displaying Parametrically Defined Surfaces, CACM, 
Vol. 23, No. 1, 23-34. 

[Boehm, 1984]. Boehm, W., G. Farm and J. Kahmann, ”A survey of curve and surface 
methods in CAGD,” Computer Aided Geometric Design 1(1), (1984), 1-60. 

[de Casteljau, 1959]. de Casteljau, P. Courbes et Surfaces a Poles, S.A. Andre Citroen, 
Paris. 

[Farin, 1980]. Farin, G., ”Bezier Polynomials over Triangles,” Technical Report 91, 
Department of Mathematics, Brunei University, Uxbridge, UK. 

[Farin, 1982]. Farin, G. ’’Smooth Interpolation to Scattered 3D Data,” Surfaces in 
Computer Aided Geometric Design, R.E. Barnhill and W. Boehm (eds.), North- 
Holland Publishing Company, 1983, 43-63. 

[Goldman, 1983]. Goldman, R.N., Subdivision Algorithms for Bezier Triangles, Computer 
Aided Design, 15, 159-166. 

[Joy, 1985]. Joy, K.I., Bezier Triangles, Technical Report CSE-85-3, Computer Science 
Division, Department of Electrical and Computer Engineering, University of 
California, Davis, California, January, 1985. 

[Lane-Carpenter, 1979]. Lane, J.M. and L. Carpenter, A Generalized Scan Line 

Algorithm for the Computer Display of Parametrically Defined Surfaces, Computer 
Graphics and Image Processing, 11, 1979, 290-297. 

[Lane-Riesenfeld, 1980]. Lane, J.M. and R.F.Riesenfeld, A Theoretical Development for 
the Computer Generation and Display of Piecewise Polynomial Surfaces, IEEE 
Transactions on Pattern Analysis and Machine Intelligence , Vol. PAMI-2, No. 1, 
January, 1980, 35-46. 

[Petersen, 1984], Petersen, C.S., Adaptive Contouring of Three-dimensional Surfaces, 
Computer Aided Geometric Design, Vol. 1, July 1984, 61-74. 


70 



